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We calculate the magnetization relaxation in the short-time regime for an 
ensemble of nanomagnets in the presence of a low frequency external AC 
biasing field at temperatures lower than the magnetic anisotropy energy of the 
individual nanomagnets. It is found that the relaxation is strongly affected by 
A C fields with amplitude larger than that of the T<i fluctuations in the nuclear 
. field. This will allow experimental probing of the nuclear spin relaxation 

£^ \ mechanism. 
^ ■ PACS numbers: 



1. INTRODUCTION 

Recent experiments! 1 ]' §i'@'§i'@ on ensembles of magnetic macro- 
molecules have shown evidence for resonant tunneling relaxation. Two of 
these experiments show evidence of tunneling in the "quantum regime", 
where only the two lowest levels of each molecule are occupied; in the "Mn- 
12" molecule,! this happens below a crossover temperature T c ~ 2K and 
in the "Fe-8" molecule,Q below T c ~ OAK. Experiments in the Fe-8 system 
have gone down to 70 mK, with no change in the relaxation characteristics 
below 0.36 K; this constitutes prima facie evidence for a quantum regime in 
Fe-8 molecular crystals. BIHI Theoretical workSBI13 on the "quantum relax- 
ation" below T c indicates that intermolecular dipole coupling and hyperfine 
coupling to the nuclear spins are necessary to explain the relaxation char- 
acteristics below T c . A number of predictions based on this theory have re- 
cently been verified experimentally;!!^' E3 this work is discussed elsewhere in 
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this volume.0 In particular, the prediction of a universal short-time "square 
root" relaxation has been verified, with a characteristic time tq which de- 
pends on both dipolar and hyperfine interactions (as well as the tunneling 
matrix element Ao). 

At present there is no direct measure of Ao in these systems, making 
it hard to verify whether the observed tq is that predicted by theory. In 
this paper we give preliminary results of a theory of quantum relaxation 
of ensembles of nanomagnets in the presence of an applied AC field. We 
show that the relaxation characteristics are strongly altered, in a way which 
should allow (i) the determination of A, and (ii) the demonstration that tq 
is controlled by nuclear spins as well as by dipolar interactions. 

We begin by considering an ensemble of nanomagnets or magnetic 
macromolecules. The high energy Hamiltonian for such systems has been 
shown to flow under reduction of temperature to a fixed point effective 
Hamiltonian given by 

H = J2 # (0) (r, {a k }) + J2 H D (ft, 7^ ) + £ Vfa, ) (1) 

r f,r k,k 

where T r is a Pauli matrix acting at molecular site r, in the Hilbert space 
of the two lowest molecular states, and a k is a Pauli matrix acting on the 
two relevant states of the k th nuclear spin; we assume that k = 1..N. The 
internuclear term V(a k , <r fe ') is usually dipolar, and \ V k k > j ~ 1 — 100kHz (and 
measurable as TrT 1 , in for example nuclear spin echo experiments). Below 
T c , Hjj is diagonal in f z (molecular flip-flop transitions are rare); it causes a 

bias £d(^*) = J2r^D(r — r)Tl at site r, which varies over a scale Ed m 0.5K 
around the sample. In the absence of nuclear spins 

#(°)(n = A < (2) 

but in general one has0 

H(°\t?, {5 k }) = A [cos{$ + £feLi a k n k . a k }f + + H.C.] 

k=l k=l 

This complex fixed point Hamiltonian contains all coupling, effects be- 
tween the molecular spin at r and the surrounding nuclear spins Of cru- 
cial importance are the diagonal hyperfine couplings u k to each a k (varying 
between ~ O.lmK — 0.5K in different systems), and the complex dimension- 
less amplitude a k for a k to flip when 7** does. l k , m k and n k are unit vectors, 
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and $ is a renormalized Kramers-Berry-Haldane phase (given by $ = ttS 

when lEifeLl \ a k\ 2 « !)• 

Typically, experiments begin by first polarizing the system, with all 
molecules aligned, and then watching the magnetization M(t) decay with 
time t. Here we assume zero applied field, but add an AC field H ac (t) = 
J2rAcos(ujt)f£ . The question that we wish to now ask is this-in the ex- 
perimentally relevant region \ Y^k=i l a fc| 2 << 1) * e - where no nuclear spins 
flip during tunneling, what will the magnetization of the crystal look like at 
short times in the presence of this AC field? 



2. THE GENERALIZED MASTER EQUATION 

We may write the magnetization of our system in the form 

M(t) =Y,J <%M(?,t,t) =Y,J ^(2P T (r,£,i) _ i) (4) 

where -Pj(r, £, t) is the normalized probability of the central spin at site r to 
be "up" (ie. in state \S Z = +S)) and in a static bias £ at time t. 

A solution for Pt(f, £, t) over timescales ~ 0(1/ oo) is rather messy. How- 
ever experimentally one is usually interested in relaxation over much longer 
timescales, in which case one can write a kinetic or "master" equation of the 
form 

Pa(r,^,t) = -W(A,u;0{P a (r,£,t) - P_ a 



-E 

r ,q 



(5) 

where P^ is the usual 2-molecule distribution function. Note that we have 
assumed that nuclear T2 fluctuations have decorrelated each pass of the local 
bias field through resonance from the others; this implies that u < NT^ 1 - 
When Amu are zero, the master equation reduces to a "static" kinetic 
equation,0 in which W(A.u:0 -> r^ 1 ^) ~ (A 2 /£o)e" l?l/fo , the nuclear 
spin driven transition rateJj'EJ We ignore higher order multimolecular terms 
P( 3 ),P( 4 ), etc. and assume approximate factorization of P^ 2 \ As before,0 
this means that the results we derive are only valid when 1 — M(t)/Mo « 1, 
where Mq is the saturation magnetization. 

Our method generalizes that given for a static applied field; we first 
solve for the rate function W(A,u;$,) for a single molecule,@ but now in an 
AC field, so that the total longitudinal field acting on 7^ is 

g ot = m+Aco S (ut) + 5&(t) (6) 
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in which 5&(t) is the rapidly varying component coming from fluctua- 
tions, and £(r) = + 6v(j) * s the slowly varying sum of dipolar fields 
and longitudinal hyperfine fields (f7v(r) = ujqM? for a system with a single 
hyperfine coupling uo and total nuclear polarization M? = J2k=i < &k > 
along the molecular easy axis). 

In what follows we assume (as noted above) that no nuclear spins flip 
during tunneling-the general results including nuclear flips will be published 
elsewhere. We also assume that the nuclear spins are in a thermal ensemble 
with kT > > ujq (it is easy to show that the AC field will drive them into such 
a high-T distribution; and all experiments so far have kT >> ujq anyway). 

In present experiments on magnetic macromolecules where the molecu- 
lar spin S w O(10), an AC field of amplitude SH = A/qubS equal to 1 G 
is equivalent to a bias amplitude A rs 50MHz ps 25mK. Experiments can 
range roug hly between 10 _5 G < SH < 100G ( depending on u). Thus we 
will assume that Af > A » A, where Af ps 0.5 — IK is the total spread in 
£ caused by dipolar and hyperfine fields; however the ratio A/£o is arbitrary. 
We then find the following results. 

(i) v4/£o >> 1 (large AC amplitude). The transition rate W(A,u;£) w 
W (C/A)Q(A 2 -f 2 ), where 

w m) = (7) 



for A- |f | » iljlA and 



Woit/A) = ^ (8) 



for A-\£\ <e /2A. 

(ii) A/£q < 1 (small AC amplitude). In this regime the dynamics are 
controlled by the nuclear T2 field; one gets, for the transition rate, 

W(A,u;0*^e-^ (9) 
so 

The essential effect of the AC field, when A >> fo> is to spread the 
resonant tunneling over a much larger energy range 2A (but at a reduced 
rate, except when |f | ps A, and the AC field itself varies slowly, so that f (t) is 
dominated by T2 fluctuations) . As the amplitude of the AC field is decreased 
the dynamics become completely dominated by nuclear T2 effects with the 
crossover occuring at A fo- 
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3. SHORT TIME DYNAMICS 

The solution for M(t) is now straightforward, following our previous 
methods JliJ For large AC amplitude one gets 

M (t) _ -m> j A _ A di ¥ —m- Wj Wo{m m 

for an ellipsoidal sample, where the internal field has Lorentzian spread 
r 0) = ^iN 1 - M (t)) and mean E(t) = cE D {\ - M(t)), and c depends on 
the specific geometry of our ellipsoid. For small AC amplitude one gets the 
square root relaxation found previously .0 

The solution to the large AC amplitude problem (10) is clearly not 
square root at short times; in fact one gets 

1 - M(t) m e- t/Tac ( 1 - M(t) « A/E D ) (11) 

where M(t) = M(t)/M ; the rate r^ 1 « A 2 /A, which is a factor E D /A 
faster than the relaxation "rate" t^ 1 ~ A 2 /Ed which one obtains for an 
ellipsoid in the square root regime .O At longer times, once the internal field 
E{t) » A, we get 



1 - M(t) « yjt/TQ ( 1 - M(t) » A/E D ) (12) 

Clearly if A/ Ed is not too small, we may never see a clear square root 
relaxation; correlations in will emerge before the square root does. 

For a non-ellipsoidal sample we find in general that the exponential 
relaxation rate will become (t*™^ " 1 )" 1 ~ A 2 /Ed = (A/Ed)^ 1 , since the 
internal fields are now spread over a large range ~ Ed in bias space. In all 
but specially-shaped samples, this will be what is observed. 

It is clear from these results that AC experiments can test the nuclear 
spin mediated tunneling mechanism, since the prediction here is that it is 
only when A > £o that significant deviations from the static square root 
relaxation will appear. This prediction is quite different from what one would 
find if we assumed that tunneling proceeded at a rate A for molecules near 
exact resonance (|^| < A); in this case the square root relaxation would break 
down for A > A. Since in the Mn-12 and Fe-8 systems, A sa 10~ 10 - 10~ 9 iT, 
this would imply a breakdown of the square root law once A exceeded roughly 
10 -8 — 10 -7 G, a very small value indeed! 

Notice also that since r^, 1 ~ A 2 /A, we have a way of determining 
the important parameter A directly in experiments, knowing A. It is clear 
that the results of AC experiments conducted near H = in magnetic 
macromolecular crystals will give a crucial test of present theory. 



G. Rose and P.C.E. Stamp 



ACKNOWLEDGEMENTS 

This research is supported by NSERC, the CIAR, and the ILL; PCES 
thanks Prof. P. Nozieres for his hospitality in Grenoble. 

REFERENCES 

1. C. Paulsen and J.G. Park, in "Quantum Tunneling of Magnetization- QTM'94" 
(ed. L. Gunther and B. Barbara), Kluwer publishing, pp. 189- 207 (1995). 

2. M. Novak and R. Scssoli, pp. 171-188, ibid. 

3. D.D. Awschalom et.ai, Science 258, 414 (1992) and refs. therein. 

4. B. Barbara et.al, J. Mag. Magn. Mat. 140-144, 1825 (1995). 

5. J.R. Friedman et.al, Phys. Rev. Lett., 76, 3830-3833 (1996). 

6. L. Thomas et.al, Nature 383, 145-147 (1996). 

7. C. Sangregorio et.al, Phys. Rev. Lett., 78, 4645 (1997). 

8. N.V. Prokof'cv and P.C.E. Stamp, J. Low Temp. Phys. 104, 143 (1996). 

9. P.C.E. Stamp, pp. 101-197 in "Quantum Tunneling in Complex Systems" , ed. 
S. Tomsovic (World Scientific, 1998). 

10. I. Tupitsyn, N.V. Prokof'ev, P.C.E. Stamp, Int. J. Mod. Phys. Bll, 2901 (1997). 

11. N.V. Prokof'ev and P.C.E. Stamp, Phys. Rev. Lett, (in press); cond-mat 
9710246. 

12. T. Ohm, C. Paulsen, preprint (Dec. 20, 1997); and see this volume. 

13. L. Thomas, B. Barbara, preprint (Jan 16, 1998); and see this volume. 

14. N.V. Prokof'ev and P.C.E. Stamp, this volume. 



